Electromagnetic surface states in structured perfect-conductor surfaces 
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Surface-bound modes in metamaterials forged by drilling periodic hole arrays in perfect-conductor 
surfaces are investigated by means of both analytical techniques and rigorous numerical solution of 
Maxwell's equations. It is shown that these metamaterials cannot be described in general by local, 
frequency-dependent permittivities and permeabilities for small periods compared to the wavelength, 
except in certain limiting cases that are discussed in detail. New related metamaterials are shown to 
exhibit exciting optical properties that are elucidated in the light of our simple analytical approach. 
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Structured metal surfaces offer a playground to yield 
remarkable optical phenomena ranging from extraordi- 
nary light transmission through subwavelength hole ar- 
rays to innovative types of surface resonances 0] ■ The 
relevant role of surface plasmons at visible and near- 
infrared frequencies was emphasized in early develop- 
ments 0,0, while subsequent studies, following pioneer- 
ing works in the field 0, l1 , demonstrated similar effects 
in plasmon-free good conductor films at microwave and 
THz frequencies |8( ■ Recently, these two regimes have 
been connected by the exciting prediction of Pendry et 
al. of surface resonances that mimic surface-plasmon 
behavior in perfect-conductor surfaces (PCS's) textured 
with subwavelength holes, and by its subsequent experi- 
mental observation 

Surface plasmons, originally predicted by Ritchie 0, 
have given birth to the rapidly growing field of plasmonics 
owing to their potential application in areas as diverse as 
biosensi ng llfi , information processing via metal-surface 
circuits |ll|. or laser technology [l^. Hence the impor- 
tance of devising new ways to achieve surface-plasmon- 
like behavior in different frequency domains (e.g., using 
phonon-polaritons in the infrared 0, Q|, or textured 
PCS's at lower frequencies 

In this Letter, we introduce a new approach to system- 
atically study surface resonances in structured PCS's. 
This allows us to provide further insight into recently 
proposed holey metamaterials , for which we find sim- 
ilar qualitative behavior and significant quantitative cor- 
rections in the surface mode dispersion relations. Fur- 
thermore, we show that this kind of materials cannot be 
represented in general by local, frequency-dependent op- 
tical constants [e(w) and except in some limiting 
cases. Finally, our results suggest a new systematics in 
the analysis of textured PCS's that is applied to related 
metamaterial designs. 

The material proposed by Pendry et al. 2] consists of 
a planar PCS perforated by infinitely-long square holes 
of side a arranged in a periodic square array of period 



d <C A and filled with homogeneous material of per- 
mittivity £h and permeability //h (^h = 1 will be used 
throughout this work), as sketched in the inset of Fig. 
^i. The assumption that the field inside the square-holes 
can be approximated by the lowest-frequency mode (i.e., 
the TEi : o mode 0), allows them to describe the bulk 
material by effective local optical constants e|| and for 
fields parallel to the surface and e± = ji± = oo for per- 
pendicular fields (PMC model Q)- In such an effective 
medium, the reflection coefficient for p-polarized incident 
light (external magnetic field parallel to the surface) is 
given by Fresnel's equation 
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where k = uj/c is the free-space light momentum, and 

/c|| and k z = ^Jk 2 — fcjj are the momentum components 

parallel and perpendicular to the surface, respectively. 
The surface resonances are signalled by the divergence of 
r p for incident evanescent waves (k\\ > k), leading to 
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with A = a 2 (the hole cross-section) and T = GA/j^/tt 6 in 
the long- wavelength limit . 

Fig. H illustrates through a representative example a 
comparison between the PMG model |2j (dashed curves) 
and a rigorous numerical solution of Maxwell's equations 
obtained by expanding the electromagnetic field in terms 
of diffracted plane waves outside the surface and square- 
waveguide modes inside the holes. Our calculated re- 
sults exhibit the same qualitative behavior as the PMG 
model in the long-wavelength limit (Fig. ^,), although, 
as expected, the inclusion of diffraction orders outside the 
material introduces new structure in the high-frequency 
region and bends the dispersion relation at the bound- 
ary of the first Brillouin zone. Closer inspection into the 
long-wavelength region reveals a sizable correction in the 
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position of the resonance towards smaller values of ku 
(Fig. QJd) , quantified in a factor-of-3 larger decay length 
of the surface state into the vacuum. The physical origin 
of this discrepancy can be understood within the analyti- 
cal approach that follows, corroborated by full numerical 
solution of Maxwell's equations. 

We consider a unit p-polarized plane wave incident 
along the xz plane (no surface modes are obtained for 
s polarization). The material will be contained in the 
z < region. In the absence of any surface structure, 
the total (incident plus reflected) field reads 



E ext (r) = -r[ik z sin(fc 2 z)x - ku cos(k z z)z] e 

k 

H oxt (r) = 2cos(fc z z)ye 
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Now, the small hole limit (a <C d, A) allows representing 
each hole by effective dipoles. This was shown by Bcthc 
[l7| for a single hole in a thin perfect-conductor screen, 
where the field scattered from the hole is equivalent to 
that generated by an electric dipole perpendicular to the 
surface plus a parallel magnetic dipole, which are in turn 
proportional to the external perpendicular electric field 
and parallel magnetic field via the polarizabilities cue and 
Qm, respectively. Parallel electric dipoles and perpendic- 
ular magnetic dipoles are forbidden by the condition that 
the parallel electric field and the perpendicular magnetic 
field vanish at a PCS. A single hole in a semi-infinite PCS 
can be described in the same fashion. In the electrostatic 
limit, Q!e ((»m) can be calculated from the electrostatic 
(magnetostatic) solution for an external electric (mag- 
netic) field, as shown in Fig. (Fig. Et>) , which involves 
only TM modes (TE modes) of the hole cavity. 

The actual self-consistent field acting on each hole in- 
cludes contributions from inter-hole dynamical interac- 
tion and symmetry considerations show that the mag- 
netic and electric dipoles must be oriented as m = my 
and p = pz. Actually, m and p depend on hole positions 
R = (x,y) just through phase factors exp(i/c||x). This 
permits writing the self-consistent relations [f| 

p = a E [E?* + G™p + G™m\ 



III 



r tt 



ZZ 

oxt i /yMM 



V 



G^m + G^p], (4) 



where GW = G™ M and G™ = -G^ E = -GS M are 
lattice sums of the electric (E) and magnetic (M) dipolc- 
dipole interactions, and i and j denote y and z Cartesian 
components. More precisely, 

Gf = E e ~ ifc "*( fc2 + ^)V 
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This inter-hole interaction is generally small for a <ti d, 
except when a diffraction order goes grazing, in which 



case the above sums can diverge giving rise to phenom- 
ena related to Wood's anomalies 0. It is near these 
divergences that surface-bound modes can exist, subject 
to the condition 
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which is obtained from the vanishing of the secular de- 
terminant of Eqs. Q). In particular, near the light-line 
in the ku — lu plane for ku > fc, one has 



Re{G^} « Re{Gyy} « Re{G^ M } 
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which comes from the divergent term of the parallel- 
momentum expansion of Eqs. JSJl (see Fig. 0a,). Further- 
more, upon inspection, one finds that Im{Gj^} = 0, and 
the remaining imaginary parts of all quantities in Eq. © 
cancel out exactly since ImjGJj'} = Imja^ 1 } = — 2fc 3 /3 
[lflj . with v = E,M. Combining these results, Eq. © 
can be approximated by Eq. with 
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Eq. ijHJ) is exact in the a -C d -C A limit, and it predicts 
the existence of surface-bound modes under the condition 
1/Rejag 1 } + 1/Rejag 1 } > 0. Calculated values of T are 
offered in Fig. |2t for various systems. The position of 
the surface mode calculated from Eq. (JSJ (see Fig. \I]p) 
differs slightly from the exact numerical result, mainly 
due to neighboring-holes multipolar interaction for a — 
0.8d (the holes occupy 64% of the surface). 

This description in terms of effective dipoles permits 
writing the specular reflection coefficient for p-polarizcd 
light as r p = I + S(m — pk\\ /k) . It is easy to see that this 
expression does not conform in general to the assumption 
of local optical constants implicit in the derivation of Eq. 
Q). In particular, for non-grazing incidence and A > (I, 
the dipole-dipole interaction can be neglected in Eqs. (@J, 
so that using explicit expressions for the fields as provided 
by Eqs. 10), and noticing the r p deviates only slightly 
from unity under these conditions, one finds 
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which is independent of fcii (i.e., of the angle of incidence) 
only if cue = 0. Otherwise, the metamaterial is non-local, 
so that the optical constants of an equivalent homoge- 
neous medium will depend on both frequency and mo- 
mentum (spatial dispersion). It should be noted that the 
neglect of cavity modes other than the lowest-frequency 
one (TEi j0 ) yields cue = (see Fig. Et), and therefore, it 
leads to an incomplete local-response description of the 
metamaterial Q. 

From the point of view of external fields, the local- 
response picture will be still maintained if |q;e| <C |q;m|) 
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so that the second term in the right hand side of Eq. © 
can be overlooked. Such metamaterials can be achieved 
by filling the holes with media of very-high | e^j 3> 1 or 
alternatively by using specific electrostatic resonances in 
the — 1 < eh < range (piling up towards eh = — 1 as 
a manifestation of their filling-material surface-plasmon 
origin). This makes cue negligible in the long- wavelength 
limit, as shown in Fig. ^ 

Another possibility consists in filling dimples rather 
than holes using moderate values of eh > 1 (e.g., £h = 10 
in Fig. |SJ) and X/d. Indeed, the high-frequency propa- 
gating modes of infinitely-long holes become resonances 
of finite width (coupling to the external continuum) in 
holes of finite depth (dimples), which are blue shifted 
with respect to the noted propagating modes owing to 
reflection at the bottom and top ends of the hole (Fabry- 
Perot picture), as illustrated by comparing Fig. [2) and 
Fig. For the choice of parameters of Fig. Ot, there 
is a resonance in om near the lowest-frequency cavity 
mode (A = 7.2a), where cke is comparatively negligible. 
This dimple resonance provides a cut-off of surface-bound 
modes 2], as illustrated in Fig. [5] Higher-energy reso- 
nances of both TM and TE nature are also observed in 
the reflectivity, giving rise to a rich structure of surface- 
bound states (Fig. EJ . 

Similar behavior could be obtained by exploiting the 
mentioned single-hole electrostatic resonances. Optical 
phonons in alkali halides yield £h < and can be 
combined with noble metals (near-perfect conductors) to 
implement these ideas in the THz domain. From Eqs. @ 
and JSJ and from the electric polarizability given in Fig. 
0J the frequency dependence of eh (a;) produces dispersion 
relations similar to Fig. El (not shown), including a cut-off 
due to the electrostatic resonances near £h = — 1. 

The surface modes resemble surface plasmons not only 
in their limited penetration into the vacuum, but also 
in the interaction that they provide between additional 
surface features like holes of larger dimensions. Indeed, 
the scattered field produced by one of such features 
in the metamaterial decays away along the surface as 
exp(i/csp-R)/v^R at large distance R. This has the same 
form as the charge distribution accompanying a surface 
plasmon launched by a localized source |2l| . in contrast 
to the exp(ifc-R)/ R far-field dependence of the interac- 
tion on unstructured PCS's 22]. Furthermore, the far 
field of a small, localized additional surface feature can 
be assimilated to the field of an effective dipole placed at 
the surface of an equivalent homogeneous material with 
the same reflectivity as the holey metamaterial. Inter- 
estingly, in contrast to the dipoles that describe the un- 
derlying hole structure, the new effective dipole can have 
parallel electric and perpendicular magnetic components. 
This introduces another handle in the design of surface 
states by using the above holey metamaterials as the base 
fabric to build metamaterials drilled by larger holes, al- 
lowing us to speculate on surface modes in fractal struc- 



FIG. 1: (a) Modulus of the specular reflection coefficient |r p | 
of a perfect-conductor surface (PCS) perforated by infinitely- 
long square holes (see inset for parameters) as a function of 
wavelength A and momentum parallel to the surface fcy. The 
surface mode predicted in Ref. |2 is shown by a dashed curve, 
(b) Detail of the reflection coefficient r p along the AB seg- 
ment of (a). 



FIG. 2: (a) Electrostatic electric-field flow lines for a hole 
drilled in a semi-infinite perfect-conductor subject to an ex- 
ternal field E cxt perpendicular to the surface, giving rise to an 
electric dipole p = aE-E cxt as seen from afar, (b) Magneto- 
static magnetic-field flow lines for the same hole subject to an 
external parallel field H cxt and leading to a magnetic dipole 
m — ctMH cxt . (c) Summary of polarizabilities for square and 
circular holes in PCS's, normalized using the hole area A. The 
circular hole in a perfect-conductor thin screen is analytical 



tures that imitate Sierpinski's carpet. 

In conclusion, we have introduced a formalism to study 
textured perfect-conductor surfaces that allows us to ob- 
tain quasi-analytical long- wavelength exact dispersion re- 
lations for surface-bound modes. We have found that 
these metamaterials cannot be assimilated in general to 
equivalent effective homogeneous media described by lo- 
cal optical constants, except in some limiting cases (e.g., 
by filling the holes with high-index-of-refraction mate- 
rial). Finally, our results pave the way towards simple 
analysis of new metamaterial designs based upon the co- 
existence of different hole sizes and distributions that can 
realize the goal of achieving tailored surface dispersion 
relations. 
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